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The conductivity tensor is introduced for the low-dimensional electron systems. Within the 
particle-in-a-box model and the diagonal response approximation, components of the conductivity 
tensor for a quasi-homogeneous ultrathin metal film and wire are calculated under the assumption 
d = Af (where d is the characteristic small dimension of the system, Af is the Fermi wavelength 
for bulk metal). We find the transmittance of ultrathin films and compare these results with 
available experimental data. The analytical estimations for the size dependence of the Fermi level 
are presented, and the oscillations of the Fermi energy in ultrathin films and wires are computed. 
Our results demonstrate the strong size and frequency dependences of the real and imaginary parts 
of the conductivity components in the infrared range. A sharp distinction of the results for Au and 
Pb is observed and explained by the difference in the relaxation time of these metals. 
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I. INTRODUCTION 

Thin-film materials, in particular, metal films are 
used widely in modern technologies including electron- 
ics. As a rule, ultrathin films are fragmentized (the is- 
land films) and consist of flat islands connected with the 
thin threads-channels [H-Q- 

Experimental techniques allow the optical characteris- 
tics in the infrared range to be measured not only for thin 
films (see, for example, but also for the specifi- 

call y g rown nanorods-antennas of the micrometer length 

In Q, the authors for the first time measured the in- 
frared conductivity of Pb ultrathin films. A decrease in 
the conductivity of the films was explained by their gran- 
ular structure. Subsequently, Tu et al. Q measured the 
optical characteristics of metal films at a temperature of 
10 K and revealed an anomalous optical transparency in 
the far-IR range. Pucci et al. [14| were the first to study 
the quantum size effects in the transmission spectra of 
lead thin films by IR spectroscopy. It should be noted 
that results of measurements have been usually inter- 
preted by experimenters in the framework of the modified 
Drude theory. A theoretical analysis of optical proper- 
ties of ultrathin films and wires is necessary, in particular, 
for the diagnostics of the nanostructure materials |19] in 
order to use them in micro- and nanoelectronics [20j. 

The important feature of the metal ID- and 2D sys- 
tems, films and wires, is an anisotropy of their electrical 
and optical properties caused by the size quantization. 
For this reason, the conductivity of the low-dimensional 
systems is represented by a tensor a a p(q, ui) which, in 
particular, determines the optical absorption. The dissi- 
pation of energy of the plane monochromatic electromag- 
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netic wave with the frequency uj and the wave vector q 
in unit volume per unit time for a nonmagnetic material 
is 

Q(q,w) = i^K^q,^) + a^(q,o;)} E a E* p , 

a, 

where E a p are the components of the electric field [2lj ]. 
However, the only value directly measurable for an ultra- 
thin film in IR range is the transmittance. 

The purpose of this work is to calculate components of 
the conductivity tensor for quasi-homogeneous ultrathin 
metal films and wires provided the condition d = Af is 
satisfied. We use the Wood and Ashcroft approach [22 1 
adapted to this case. The main advancement is the pro- 
cedure of the accurate determination of the Fermi level 
for a film and a wire of such a thickness with taking into 
account the size oscillations. The transmittance of the 
ultrathin films is also calculated in order to compare the- 
oretical results with experimental data [3, [j| . 

To the best of our knowledge, the oscillatory behavior 
of the Fermi energy in confined (2D) electron gas was 
studied for the first time by Sandomirskii [23| . Later 
calculations were performed for thin films, spheres and 
wires [241 - I341 ]. In the present paper, in order to describe 
the Fermi energy behavior in low-dimensional metallic 
systems we use an elementary one-particle analytical ap- 
proach [32| . 



II. CONDUCTIVITY TENSOR 

A film of thickness L (or a wire of radius po) compara- 
ble in magnitude to the Fermi wavelength of an electron 
in an infinite metal (A F ~ 0.5 nm) will be referred to as 
the ultrathin film or wire (see Fig. 1). The longitudinal 
sizes of the sample are assumed to be considerably larger 
than the film thickness: L« fl,4 (or po "C C for wire), 
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which leads to the pronounced quantization of the trans- 
verse component of the electron momentum. This re- 
sults in the formation of subbands, i.e., groups of energy 
levels corresponding to the same value of the transverse 
momentum component. 



A response of an electron gas to the electromagnetic 
field E = E cxp[i (qr — Lot)] may be determined in a 
linear approximation by the density matrix technique. 

For the induced current one can obtain 1221: 
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Figure 1: Choice of coordinates. 

where io = \i), \j) are the wave functions of the 

initial and final electron states corresponding to energies 
£i and £j \ Sij — £, — Sj] fj and fj are occupation factors; 
Q is the volume of sample, m e is the electron mass, — e 
is the electron charge, p is the momentum operator. 

We describe conductivity electrons in metal films 
and wires within the framework of the particle-in-a-box 
model. In many cases, this model turns out to be quite 
productive for metals with high conductivity (35j . In the 
case of the low-dimensional electron systems, a potential 
box model includes the key feature of these systems - 
a confinement of electrons inside a region with certain 
small dimension. The model remains applicable up to a 



certain critical size when conductivity vanishes. As for 
the depth of a box, this parameter remains important as 
long as electron emitting is permitted. 

In general, the model of electron gas in a rectangular 
potential box is a good initial approximation, while var- 
ious details (particularities of structure, impurities etc.) 
can be taken into account by introducing corresponding 
corrections. 

An electron system in low-dimensional structures is 
anisotropic, and its characteristics can be expressed in 
tensor form. The tensor origin of the conductivity be- 
comes obvious after converting the expression (fTJ) into 
the form 

ja (k, w) = ^2 a a{s (k, q, u)Ep (q, u), 

where a, /3 = x,y,z and a a p is the conductivity tensor. 

It is not difficult to demonstrate that the conduc- 
tivity tensor is proportional to 5 Q ^(5k,q, where 5 a p or 
<5k, q = {1, k = q; 0, k ^ q} is Kronecher's symbol, 
for macroscopic samples with the wave functions of the 
kind SI -1 / 2 exp (— pr/fi). This implies that all the Fourier 
components of the current, except the one with k = q, 
are equal to zero. Of course, this is not the case for ul- 
trathin films and wires, but the component with k = q 
is still dominating. At the first step, known as the di- 
agonal response approximation, this component only is 
taken into account. We then find 



a^(q,q, u) = ^S a , + £ AJl ( (j\ e~^p a \i) - hq a (j\ e"W 

m e uj\l mjujil ^— 4 — Hu) \ 2 J 

x U\ \ j} + 1 % (i | e w ^ = a^cj). (2) 
I 



Here the relation ^\ fi = Nis used with N equal to the number of the conductivity electrons. 



3 



Over infrared region, the condition qL, qpo <§C 1 is sat- 
isfied allowing us to express the conductivity tensor in 
terms of the according small value. 



III. FILM 

It is assumed that the conduction electrons of the film 
are located in a rectangular potential box V(r) with a 
depth Uq < 0, so that the box shape reproduces the film 
shape (see Fig. 1), and 



\U \=4 + W , e F = ^{^ 2 n) 2 '\ 
2m 



(3) 



Here Wq, e F , and n are the electron work function, the 
Fermi energy and the electron concentration for a bulk 
metal, respectively. 

The unperturbed states of the film are described by 
the wave functions 

*mn P (x, y, z) = -l=i; m {x)e 2 ^y' a e 2 ^' z l\ (4) 
V ab 

where n, p = ±1, ±2, . . . and m = +1, +2, . . .. The sub- 
script m numbers the subbands. The wave functions 
ip m {x) are represented in the following form: 
for even values of m, 



( C m sink xm x, -L/2 < x < L/2, 
1> m (x) = { (-lY m ^ +1 B m e-^ x , x > L/2, 
(-l) m / 2 B m e K ™ x , x < -L/2, 



and for odd values of m, 



C m cosk xm x, —L/2<x<L/2, 

(•_ 1 )(m-l)/2 jBm g- Km X j x > L /2 : 

(-l)( m - 1 )/ 2 B m e"» iB , x < -L/2, 



(5) 



(6) 



2/sv 



B 7n — C, 



k 



xm c K m L/2 



Z + K rn L Kq 

Here, C m is the normalization factor, k xm are the roots 
of the equation 

k xm L = -2 &rcsm(k xm /k ) + nm, (7) 



where n m = \/k 2 — k xm and hko = y/2m e \Uo\ (see Ref. 
32]). 

In this section, we focus on optical transitions between 
subbands accompanied by changing the transverse com- 
ponent of the electron wave vector k xm . These transitions 
contribute to the a xx component of the conductivity ten- 
sor. Since qL <C 1, we have in zero approximation 

i e 2 



x JV + 



J_yji-fj 



01 



>-io(q y y+qzz') 



Vx\l) 



(8) 

Dividing by — fko in the sum and interchanging i and 
j for the second term appeared after this dividing, ex- 



pression ([8]) can be transformed into 

•,2 



0\j 



01 



N+—T] J% 2 



Since 



0> 



-io{q v v+qzz)f. 



Px\i) 



(9) 



Px\i) = (™>'\Px\m)S qy ,k yn -k yn ,8qz,k zp -k zpl , 



and in view of the fact that \k xm — k xm '\ 3> q, further 
simplifications are possible: 



( 



N+—Y] _^ £ —' 



■Km'lp^lm)! 2 



V 



m,rn 
n,p 



(10) 



Here the occupation factor is approximated by the step 
function f mnp — 9(ef — £ m np), where £f is the Fermi 
energy for nanofilm, e mm ' = h 2 (k 2 m — k 2 ,)/2m e . Using 
Thomas- Reiche-Kuhn sum rule (see Ref. [22j), we rewrite 
(Ha as 



2i e 2 h 2 uj ^ fmn P \(m'\p x \m)\ 2 



(11) 



and then one can obtain corresponding component of the 
dielectric tensor 



= 1 



47TZQ 



(12) 



The matrix elements of the momentum projection op- 
erator p x — ihd/dx from (QJ - © are 



|(m'|^|m)| 2 = {l-(-ir +m '} 



8^ k xm k xm ,K m K m i 



{k 2 xm ,-ki m y{2 + K m L){2 + K m ,Ly 



(13) 



The broadening is introduced in a manner proposed 
by Mermin (3|| . As a result of this procedure the tensor 
components a xx and e xx get both real and imaginary 
parts: 



Re 



Im g 



\LJ tt \h2a ) 

4\3 ao fc2 /l e 2 \ 
LJ 7T \h2an) 



(+)' 



(14) 
(15) 



Ree xx 



1 - ( z) 4 ^-» Ime * 



4x4 W 

l) ^k! H ^ 



(16) 
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where 

?X r ^ x , T 2 k k ,k 2 k 2 (k 2 -k 2 \Uk 2 - h 2 \ 2 □= A- 4 in vH 

^"i^i^ 1 1 j 1 (2 + Km I)(2 + Km ,L)(fcL, - fcL) 3 {{(^L " ^J 2 " ^ + 7 4 } 2 + 4fc*7 4 } ' 1 J 



Here 7 = ^2m e /hT, r is the relaxation time, hk u 
\/2m e hio, a,Q is the Bohr radius and 



— arcsm 

7T 



(-) 







~Lk ~ 










IT 



■ 1. 

(18) 

Square brackets in (TTS)) and in the text below indicate the 
integer number. Instead of the summation over n and p 
in (fTTj) . (fT"2]) we perform integration. 

In order to use Eqs. (fl4)) - (Tl8)) in calculations, it is 
necessary to supplement them by the relation determin- 
ing the Fermi energy of film [32j 



k'^ 



1 / ™ F \ 

= 2nnL + *4J 



(19) 



The relation (fig]) together with Eqs. and CED de- 
scribes the size-dependent Fermi level in ultrathin films. 

In the case of a film, transmittance is a quantity, which 
is directly measurable 



TR = 7/7 , 



(20) 



where Iq and I are intensities of a wave at surfaces x = 
—L/2 and x — L/2, respectively. 

For a film of thickness L the transmittance may be 
estimated as: 



corresponds to the longitudinal motion of an electron. 
The subscript p numbers values of z— component of its 
wave vector. The angle part of the wave function 



® m (¥>) 



1 



/2vr 



has to satisfy the periodicity condition 
$ m (if + 2tt) = $ m (tp) 



(25) 



(26) 



from which follows the eigenvalues spectrum m = 0, ±1, 
±2,.... 

The radial dependence of the wave function is de- 
scribed by the Bessel functions of an integer order 



where 



Rmn (p) C mn I m (fcmnP) j 

V2 



PO \I' m (k mn p Q )\ ' 



(27) 



(28) 



Here k mn = a m n/po, where a m „ are positive roots of the 
Bessel function of the m-th order 7 m (£), n = 1,2,.... 
The prime marks a derivative with respect to £. 

In the next Section, for the case of a wire, we obtain 
the relation similar to Eq. (TTJ 



TR = exp{-r](oj,L)L}, 



(21) 



where the absorption coefficient rj should be calculated 
by using (JTBJ) and the formula 



7] = — Im^e(w, L). 



IV. WIRE 



(22) 



The simplest model for an ultrathin wire (see Fig. 1) 
is to consider it as a cylindrical potential well V(p, z) 
of infinite depth. The length of the well £ is assumed 
to be much larger than its radius po- The conductivity 
electrons are described by the wave functions of the kind 

Iprnnp (p, ¥>, Z) = Rrnn (p) $m (v) ^p (?) • (23) 

The function 

1 



Z p (z) = 



(24) 



A. The Fermi energy 

We start from the expression for the energy of an elec- 
tron 



2m, 



(* 



2 

mn 



(29) 



where k mn and k zp are the eigenvalues of transverse and 
longitudinal components of the electron wave vector, re- 
spectively [20(1 ■ The electron states in a wire correspond 
to points (k mn ,k zp ) on the k± k zp half-plane (k± > 0). 
Since the spectrum k zp is quasicontinuous (£ S> po), 
these points form a system of straight lines k± — k mn . 
The occupied states distribute on intercepts cut off by 
the semicircle of radius fcp (see Fig. 2). Density of the 
electron states at the intercepts is equal to C/tt. 

The total number of the occupied states (equal to the 
number of the conduction electrons in a wire) is 



N = 2- 



jU2 _ 1,2 
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Figure 2: Geometrical diagram of electron-state filling in 
quantum wire. 



Taking into account that TV = nf2, we then obtain the 
equation for the computation of the Fermi level fcp in an 
ultrathin wire 



' u m . n. 



U2 



(30) 



The electron concentration n is assumed to be the same 
in a wire and in a bulk metal. The summation should 
be performed over all numbers to and n satisfying the 
condition 



< k ¥ 



(31) 



The size dependence of the Fermi level in ultrathin 
films and wires has an "oscillatory" form. In order to 
determine a magnitude of the variations, let us evaluate 
the averaged (smoothed) size dependence. In this case, 
an averaging means a replacement of summation in Eqs. 
(QH and ([201) by integration. 

For a film, we use the Euler-MacLaurin summation for- 
mula (37l ] . in which it is enough to take into account the 
first two terms. Allowing m F to take any value (not only 
integer, m F ~ Lk-p/n) and neglecting the corrections for 
a finite depth of a potential box, we obtain 



fc F /fc F «l + 7r/(4fc F L), 



(32) 



where k F is the Fermi wave number for a bulk metal. 

In the case of a wire, it is hard to estimate directly the 
size dependence of the Fermi level because it is impossible 
to express explicitly the roots of the Bessel functions. 
However, the averaged size dependence k F (po) can be 
obtained in an indirect way. 

Let us rewrite (f3"Ul) as 



(k F p ) 3 =6j2\f(k F Po) 2 -al m . 

m,n 

Here fcp = (Stt 2 ?^) 1 ^ 3 , and the relation 



(33) 



was used. a mn are positive roots of the Bessel function 
of order to = 0, ±1, ±2, . . . and n = 1, 2, Assuming 



the a mn = a(m, n) function to be continuous, we turn to 
the integration: 



(fc^o) 3 = 12 yy ^Thpj) 



a 2 (m,n) dmdn (34) 



(to > now). Limits of the integration are determined 
by the condition a(m,n) < k F po. 

The left-hand side of Eq. (j34l) tends to zero with po —> 
0. At the same time, the right-hand side tends to zero 
only if k F p — > a i. Hence, the averaged size dependence 
is 

k F (p ) sa a 01 /p 

for the small values po- For the large values po, the in- 
tegrated expression is k F po and the area of the region 
of integration J J dmdn is proportional to k F p\. Then, 
comparing the left-hand and right-hand sides, we find 
that fcp(po) ~~ ► const with p — > 00. Accepting the con- 
stant to be fc F , we finally obtain 

k F /k F = l + a 01 /(k F p ), 

where aoi ~ 2.4048. 

In the next subsection, to calculate conductivity com- 
ponents, we use the size-dependent Fermi energy e F 
found from the exact expression ((30|) . 



B. Components of the conductivity tensor 

Let us consider the case when a wave is directed nor- 
mally to the axis of a wire (see Fig. 1). The wave vector 
is then located in the x—y plane, i.e., q z — 0. Orientating 
the x— axis along the wave propagation, we get q y = 0, 
qr = q x x ~ po/X <C 1, and e ±,oqr »1± ioQxX. 

In zero order of the expansion a a fs in terms of po/X, 
the expression ^ takes the form 



i e 2 N i e [ 

&a/3 = 7^ 6 aB-\ 5— 

TO e WS2 TO„Co 



\ " fi fj 



01 Pa \i) {i\pp \j) 



(35) 

Following the procedure, which led us to Eq. ©, we have 
i a e 2 N 



&af3 



fi 



! 01 Pa \i) {i\Pp \j) + 01 Pa (i\PP U) j ( 36 ) 



Using f|24[) - fj28[) . after rather cumbersome transfor- 
mations (see Appendix A) , the matrix elements of various 
projections of the momentum operator can be written as 



01 Pa 



— lo^A ,h fl G, ^ 
~2 $pp' kmnC'mnG Y-f-) ; 



(37) 
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J { _) T S m +i,m>J(+), (38) For other diagonal components, the expression ([36]) can 

be easily transformed into 



po 



Because of a specific form of (j\ p z \i) the sum in Eq. ([3"5 
becomes zero if a = z or j3 — z. Hence, 



^xz.zx.yz^zy — 0; ® zz 



ioe 2 n 



ioe 2 n 2ioe z 



E 



fi &i 



m, e uj miujQ, s% — h 2 cu 2 



I01«r, (40) 



where the subscript a = x,y. After a substitution of the 
(39) matrix elements ([5T1) into P0|) . we find 



J xx,yy 



ioe 2 n 
m e uj 



m e uj 



, fmnpk mn C t 



(k 



2 

ran 



(") 



q / j J rnnp^rnn^mn } / 7,2 
52 m.n 



+ 



(fc 



2 

ran 



k 2 



( + ) 



^2 _ 1,4 1 (1.2 _ 1.2 \ 
l,ri) ""w I'-mti "-m+ljTi'/ 



fc 4 , 



(41) 



r 



where 



f„ 



1 1.2 

x ' "'mn "T" ^zp ^ ^F' 
n t 2 , 7,2 > 7,2 



An expression for the non-diagonal components a xy 
and a yx follows from (|36[) 



m 2 ujfl ' ' 



(42) 



The upper sign corresponds to a = x, /3 = y, and the 
lower one to a = y, /3 = x. 

The axis symmetry of the problem is reflected by the 
fact that in (j35p and (|3G[) the summation is performed 
over positive to and ml as well as negative ones but the 
same in absolute value. The analysis of the expressions 
([3"T]) and ([3"5|) based on the properties of the Bcsscl func- 
tions [38j 

^(— m)n kmni -^—rriis) ( 1) -^m(0 

reveals a different behavior of the matrix elements when 
changing together to — > —to and m' — > —to', 

01 P* |») -> - (j\Px \i) , (i| p y \i) -> 01 A/ I*) • ( 43 ) 

This causes the terms in (|42|) to cancel pairwise, and we 
then find 



' xy — u yx 



0. 



(44) 



Thus, all non-diagonal components of the conductiv- 
ity tensor vanish in zero approximation of the expansion 
in terms of po/X. However, in linear approximation the 
result is different. We take account that terms, which 



contain lead to the vanishing of the sum. Therefore, 
in this approximation, components a zx and a zy have a 
form 



miujil 



5> 



1.3 



! {j\xp z \i) (i\pp \j) + {j\xp z \i)* (i\pp 



Eij — hiij 



Sij + has 



where ft = x,y and for the matrix elements see Appendix. 

An analysis, similar to the one, which resulted in Eq. 
(1431) . gives 

01 x Pz I*) -> - (j\xp z \i) , 

01 x p* K) (ilPv li> -> - 01 ^ I*) b') > 

01 x Pz \i) (i\Px \j) -> 01 ^P* |t) (i\p x \j) . (46) 

Hence, to linear order in po/X we have a zy = but a zx ^ 
0. Using Eq. (|38]l. the relation 



01 x ^ I*)* (*|Pai b')* = - 01 x ^ I*) (i\Px \j) 
and Eq. (jA3[) . we derive 
2i q x e 2 



Ml 



^2 franpkzpC 2 mn {T(^ - J r (+)), (47) 



m.p 



where 



I 'a 



p) I m (k mn p) p 2 dp 



- 



(Ul _ 1.2 \2 _ L.A 



Dissipation is introduced by the substitution w — > 
u> + iq/t in expression for conductivity. When r = 0, 
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the diagonal components of conductivity are imaginary. 
Since the remaining components of the tensor vanish in 
zero approximation, dissipation is also absent {Q = 0). 
In general, dissipation is small for optical frequencies in 
which we are interested (cj 3> 1 /t) . 

Substituting u — > to + io/r in (|39[) . after straightfor- 
ward transformations we obtain the Drude formula 139] 



(which may be obtained from (1491) and (|50p under con- 
dition that the peaks are well separated). Taking into 
account that 



E 1 



2C 



h 2 — h 2 
fop ft m 



2C 



(52) 



<J zz (uj) = cr(0) 



1 + IqLOT 
1+W 2 T 2 ' 



(48) 



where <j(0) = e 2 nr/m e is the static conductivity. Thus, 
the component u zz {uj) is associated with the classical 
conductivity. Other diagonal components (|4T)|) can be 
represented as 



<J aa = <J ZZ {1 + S(cj,p ,£)}, 



(49) 



where 



_ 2 ^ fay (4 ~ ^ 2 + 2ft 2 ^o/r) 2 
iVm e 4f (4 - ft 2 ^ 2 ) 2 + 4^ w 2 /r 2 IWaWI 

I, J J 

(50) 

and a — x,y. 

After interchanging subscripts i and j, terms of the 
sum ([50| reverse their sign. As a result, 



and 



S 



(e 2 



4/i 4 w 2 /r 2 



7Vm e ^— ^ (e 

£i<EF 

e 3 ->£f 



e« (4 ~ ^ + 2ft2 ^o/r) 2 



ft 2 w 2 ) 2 +4^ 4 o; 2 /t 2 



(51) 

Here < 0, i.e. only transitions coupled with absorp- 
tion participate in the conductivity. It is important to 
remark that ImS < for any frequency. Since in the 
optical region the real part of the component a zz can be 
ignored and its imaginary part is positive, it follows from 
(|49"]) that Re a xx> yy > and Q > over the whole region. 

Let us compare in magnitude components of the con- 
ductivity tensor. For Au, the frequency Hlo = 1 eV, dis- 
sipation h/r = 0.02 eV we find <r(0) = 4.6 x 10 17 s" 1 , 



a(0)/u 



10 
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We use below the value 



e 2 /2aoh = 2.0 x 10 16 s 1 as a unit of the conductivity. 
Then \a zz \ 0.5. 

We can now estimate, for example, height of peaks in 
Kea xx . We use relationships 



and 



ImS ps - 



Reov 



KNm, 



|cr Z2 |Im5 



|(m + l,n'\p x \mn)f^2 1 



\(m + l,n'\p x \mri)\ 2 oc k 2 mn = ^h 2 k% 2 (53) 



and, using (j3]), we have 



ImS 1 = 



3hr 



For r = 2.1 x 10~ 14 s" 1 (Au), d = 2p = 2 nm, we find 
ImS* = — 1, Reo^ = 1. In macroscopic limit po — ^ °o we 
find that Re a xx — and Im a xx = Im a zz , as we have 
expected. 

Comparing (f36|) with (|45[) . one can obtain |cr za; /o'a; a: | = 
g^po- F° r A = 10 3 nm, d = 2 nm we have \o zx ja xx \ = 

io- 2 . 



V. RESULTS AND DISCUSSION 

The difference between our approach and the theory 
[22I ] is associated with peculiarities of a electronic levels 
distribution in films/ wires of nanometers thickness, when 
the d = Af condition is satisfied. 

In this number of subbands, formed as a result 

of the size quantization, is small, while the contribution 
of each of them in the sum (|17|) is significant. Opposite 
assumptions are made in [22| : a characteristic size d is so 
large that the number of subbands is much larger than 
1. Then, the separation between neighboring subbands 
(with numbers m and m+1) is small, while contributions 
of individual items in the sum [22| are not significant 
anymore, so the summation can be substituted by inte- 
gration, as usually done in the case of a quasi-continuous 
distribution. Discreteness, coming from the size quan- 
tization, manifests itself only weakly. The Fermi level 
in films and wires with small thickness noticeably differs 
from the Fermi level of a bulk metal (30% difference for 
a wire of 1 nm diameter, see Fig. 3). In (2^, the authors 
use the Fermi level of a bulk metal to find a number of 
subbands, while we take into account the size dependence 
of the Fermi level, when determining a number of sub- 
bands. For a few nanometers thickness, these numbers 
are found to be different and, because they are small, a 
noticeable divergence in results is revealed. 

Finally, theory [22j was developed in order to apply it 
for an isotropic composite medium. Therefore, from the 
very beginning, the direction of an applied field was con- 
sidered as a preferred one. In our approach, anisotropy of 
the metallic ID and 2D systems (wires and films) is taken 
into consideration, their conductivity and dielectric func- 
tion are assumed to be tensors that allows a response of 
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Figure 3: Reduced size dependences of the Fermi energy of 
wires and films vs diameter d = 2po and thickness L, respec- 
tively. 



wires and films to be determined for any orientation in 
an external field. 

Anisotropy as well as discreteness manifests itself much 
stronger in systems with the small characteristic size d = 
Af- Under the condition d 3> Af, both our theory and 
theory [22| lead to the same results. 



A. The Fermi energy 

Fig. 3 demonstrates the size dependence of the Fermi 
energy for films and wires of Au and Al computed from 
Eqs. (IT^|) and (|3Tj|) . The size dependences have an "os- 
cillatory" form. In contrast to the Fermi energy of a 
film [32[ , the size variation of the Fermi energy of a wire 
seems to be random. Input parameters for calculations 
were taken from Ref. [32 1 . 

In the case of a film, cusps on the size dependence 
(i.e., the jumps of the derivative dep/dL) are distributed 
nearly regular with the approximately constant period 
AL « ir/k F . The cusp on the size dependence of a wire 
appears each time when the increasing radius po reaches 
the value Poim'n') f° r which the condition pip is satisfied 
by one more pair (ml ,n')\ 

a m'n' = ^FPO(m'n') ■ 

Distance between the neighboring cusps is 

Ad ~ 2 (a m 'n' CLmn) /k F , 



where a mn is the root of the Bessel function closest to 
a m i n r in value. Roots of the Bessel functions of different 
orders mix up so that Ad varies, at first sight, randomly 
with change in size. 

The oscillations of the Fermi energy in a wire of diame- 
ter d and in a film of thickness L are similar in magnitude 
if d = L. As in the case of a film, the " period" Ad and the 
amplitude of the oscillations tend to zero with increasing 
diameter. 

Characteristic properties of the size dependence of the 
Fermi energy for various metal wires (and various metal 
films too) may be explained exclusively by different value 
k F . As compared to the Au wire, for the Al wire, the 
scale Ad of the oscillations is finer, the amplitude of the 
oscillations and the averaged value ep/e F are smaller. 



B. Film 

The specific feature of the optical characteristics of 
thin films is the presence of peaks associated with the 
optical transitions between the subbands. The size effect 
manifests itself in a change of the number of peaks, their 
position, and the spacing between them. 

The positions of the peaks is determined by the ap- 
proximate expression tkv mm ' ~ hojo\m' 2 — m 2 \, where 
m and m' are the numbers of subbands between which 
the transition occurs and huiQ = n 2 h 2 / (2m e L 2 ) = 0.34 
[cY]/L 2 [nm 2 ]. The frequency range under consideration 
lies in the infrared and visible spectral ranges. The lower 
limit of the range (fiwi2) corresponds to the beginning of 
the optical transitions between the subbands. The upper 
limit of the frequency range is the electron work function 
W of the film. The estimates can be made with the work 
function Wq for infinite metals Au and Ag. 

The calculated real and imaginary parts of the conduc- 
tivity component a xx for the Au films 2 and 6 nm thick 
are presented in Fig. 4. For ultrathin films, the number 
of subbands completely or partially occupied by electrons 
is small: nip s» 2L/X F . Therefore, the number of peaks 
is small as well. For the film of the thickness L = 2 nm, 
the peak at fiwi2 ~ 0.25 eV corresponding to the lower 
limit of the frequency range is clearly seen. The peaks 
that represent the transitions between the neighboring 
subbands with the numbers m and m! = m + 1 are lo- 
cated to the left of the maximum height peak observed 
at the frequency hoj ma _ x ~ hcuQ(2mp + 1). This frequency 
corresponds to the transition between the subbands with 
the numbers m = mp and mf = mp + 1. The spacing be- 
tween any two neighboring peaks is identical and approx- 
imately equal to 2hw c . As the film thickness L increases, 
all peaks shift toward the left, the spacing between peaks 
decreases, and they begin to merge together. 

The overlapping of the peaks becomes significant when 
the spacing between them is equal to their width. The 
peak width is determined by the dissipation mechanisms 
and is approximately equal to 2H/t. Peaks for the film of 
the thickness L = 2 nm are clearly distinguishable (see 
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Figure 4: Frequency dependences of the real and imaginary 
parts of the film conductivity component (in e 2 /2aoh units) 
calculated by Eqs. (|f4[) and (|f 5|) (solid lines). Dashed curves 
correspond to the results of calculations by formulas (68) from 
work 12211. 



Fig. 4), but for the thickness L — 6 nm the peaks disap- 
pear completely. (It should be noted that the results of 
our calculations appear to be weakly sensitive to a change 
in the relaxation time r within one order of magnitude.) 

It can be seen from Fig. 4 that, as the film thickness 
decreases, the discrepancy between the results of calcu- 

from Ref. 



lations by Eqs. Q3J), and by (68) from Ref. [22 1 
increases and becomes substantial. This discrepancy is 
associated with the fact that relationships (|T4l and (fl~5|) 
were derived with allowance made for the dependence of 
the Fermi energy on the film thickness kp(L) and the ex- 
act calculation of the number mp of occupied subbands. 
In Ref. [22j, the number mp was calculated by the pro- 
cedure which gives an error ±1 for films with thickness 
L ~ Ap . This is an essential error because the number of 
occupied subbands for such thickness is small and ranges 
from 2 to 6. 

The calculated frequency dependences of the transmit- 
tance for Au and Ag thin films of different thicknesses are 
compared with the experimental data in Fig. 5. When 
comparing the results of our calculations with experimen- 
tal data, it is necessary to take into consideration that 
our definition of the transmittance is different from that 
which is normally used by experimentalists. 

In the transmittance (|20l) the value / is the same in 
both cases: This is an intensity of radiation, which comes 
out from the film through the surface x — L/2. Experi- 
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Figure 5: Frequency dependences of the film transmittance 
calculated by Eqs. (I14|l - (|22|l using e xx (solid lines) and e zz 
(dashed lines, left hand scale). Opened triangles [5|, circles 
[J], and square [9( indicate the experimental data (right hand 
scale) for the Au and Ag films. Solid triangles and circles 
indicate the recalculated experimental data (left hand scale) 
in according to Eq. ((54}. 



mentalists take Iq as intensity of radiation incident onto 
the surface of the film x — — L/2. Of course, certain 
fraction of the radiation penetrates inside, while the re- 
maining part is reflected. We don't consider reflection 
and assume Iq to be intensity of radiation, which comes 
into the film through the surface x — — L/2. To make a 
comparison with the theory, experimental values of the 
transmittance TR oxp are recalculated by using the for- 
mula 



TR = TRoxp 



(54) 



1-R' 

where R is the value of reflection coefficient obtained by 
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parts of the film dielectric function calculated by Eqs. (|16[) . 



measuring under the same conditions as TR cxp 0, @. 
The results of the recalculation are also presented in Fig. 
5. 

The absorbance ry is determined by the functions 
Imf(fiw) and Ree(fiw) according to expression (j2"2")l . The 
frequency dependences of these functions exhibit a dif- 
ferent behavior (Fig. 6). 

Unlike the function Im e(Hui), the function Re e(haj) has 
not only pronounced resonance maxima but also minima 
shown as inverted peaks. The height of both peaks in- 
creases with an increase in the frequency, so that, eventu- 
ally, one of the inverted peaks intersects the abscissa axis, 
and the function Re e(fiw) becomes negative (in contrast 
to the function Ime(fiw) that is always positive in sign). 
The minimum transmittance should be identified with 
the minimum of the function Ree(huj), which is located 
in the vicinity of the frequency hio mstx — hLj^{2m-p + 1). 

At frequencies Huj > fiw max , the absorbance is deter- 
mined only by the real part of the dielectric function: 
i] w (2w/c)-\/|Ree|. It is easy to check that the ab- 
sorbance tends to a specific constant value with an in- 
crease in the frequency. The transmittance TR (|2"Tj) 
is characterized by the same tendency. This tendency 
can be clearly seen in Fig. 5. The peaks associated 
with the transitions between far subbands are clearly dis- 
tinguished against the background of the monotonic in- 
crease in the transmittance. In particular, the transition 
m-p — 3 — > m-p manifests itself at fku ss 1.8 eV (L = 4 
nm). 



To the zero order in L/X, diagonal components of the 
dielectric tensor only are not equal to zero. Solid lines in 
Fig. 5 represent transmittance computed by Eqs. (j2"Tj) 
and (|22p using e xx . The transmittance indicates a change 
in a normal to surface component of an electric field of 
a wave passed through a film. It is this component that 
causes optical transitions between subbands formed by 
the size quantization. Dashed lines represent transmit- 



tance calculated with using e zz (e 2 



j) which shows a 



weakening of a parallel to surface component of an elec- 
tric field. Such a transmittance is observed at normal 
incidence of radiation onto a film. 

In the region ftw > ftw ma x, a mechanism of dissipation 
(i.e. value of the relaxation time r) affects the trans- 
mittance weakly. Dissipation manifests itself noticeably 
only in vicinity of the minimum of transmittance. Thus, 
discrepancies between theory and experiment at these 
frequencies could not be explained by either a mecha- 
nism of dissipation or an orientation of a film in field. A 
remarkable feature, as seen from Fig. 5, is a noticeable 
exceeding of measured transmittance over computed one 
(except of the last section, where there is a good agree- 
ment). This implies that the discrepancy between the 
theory and experiment can be attributed to a large non- 
homogeneity in thickness and especially an absence of 
continuity, i.e. it can be explained by the presence of 
regions of a substrate without coating. 

Let's name the ratio of coating area to substrate area 
by the coating coefficient p and denote coating thickness 
as L' . The mean thickness which is usually determined 
by experimentalists from the mass of a film and substrate 
area is L = pL' . Transmittance of a "holey" film of 
thickness L is 

TR(L) = 1 - p + pTR(L') = l-p + pTR(L/p). 

Reducing p, it is possible to increase transmittance up 
to 1. For example, transmittance of a film with thickness 
L = 4 nm and coating coefficient p ~ 0.45 at the fre- 
quency ftw — 1 eV is equal to 0.75, i.e. discrepancy with 
the experimental value is twice lowered (see Fig. 5). 

In the frequency region hw < 0.5 eV, there are experi- 
mental data for transmittance of thin films of Au Q and 
Pb d, 03] . Unfortunately data for reflection are absent 
hence the recalculation of experimental results like that 
represented above is impossible. However, taking into 
account that this recalculation leads to an increase of 
transmittance value, we guess experimental data for Au 
@ to be in agreement with our calculations under as- 
sumption of a normal incidence of radiation onto a film 
(Fig. 5). Moreover, the recalculation can change a type 
of dependence on frequency for transmittance, i.e. a ris- 
ing can be replaced by a falling after the recalculation 
(see Fig. 5 for Au, L = 5.8 nm). This could explain 
why within frequency interval (0.2, 0.5) eV dependence 
of transmittance on frequency is decreasing or absent at 
all according to our calculations whereas an increasing 
is observed in experiment [14} . As to value of transmit- 
tance for Pb films, it is difficult to compare calculations 
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with experimental data because results ol different ex- 
periments vary essentially. Thus transmittance value for 
Pb film of 4 nm thickness is given in [9] as 0.79 at 0.05 
eV but it is 0.12 only at the close frequency 0.12 eV fol- 
lowing 14] (in addition, transmittance diverges in value 
1.5 times for different technologies of film coating [5). 



C. Wire 

The frequency dependences of Re a xx and Im a xx for 
the Au wire of diameter 1.6 nm are presented in Fig. 
7. For such a small diameter, the peaks corresponding 
to the transitions between levels of the size quantization 
(subbands) manifest itself clearly. In spite of the rather 
complete spectrum k mn , position of the peaks is well pre- 
dicted. 

Let us find, for example, the position of the peak in 
Recr^a; which has the maximum height. The height of the 
peaks is proportional to | (m + 1, n'\p. x \ mn) \ 2 1- The 
matrix element has the maximum magnitude at nf = n 
because under this condition the integral in (|38[) takes on 
the maximum value. Further, by using (|52[) and (|53[) , it 
is easy to determine that the maximum height is realized 
at to = 0, n' = n = ?if. For the diameter d = 1.6 nm, 
the number n-p w k^po/ir is equal to 3. As a result we 
have 



te (ki 



^0,3) 



h 2 (al 



2m„ 



2TO e/ og 



1.40 eV. 



This is in good agreement with numerical calculations 
presented in Fig. 7. The height of the peaks in Fig. 7 
also confirms our estimation. 

Fig. 7 demonstrates the important fact Re a xx > 
over all frequency range. In contrast to this, Im a xx is a 
variable in sign function of the frequency. 

The frequency dependences of Recr 2a ; and hna zx are 
also presented in Fig. 7. As it was expected, the position 
of the peaks is identical both for Re a zx and Im a zx but 
the height is one order smaller in the first case. 

Comparing the upper and lower parts of Fig. 7, we 
can trace the size dependence of the conductivity for ul- 
trathin metal wires. When d increases, the peaks shift 
to the left with displacement equal to Aw = u/ — to = 
uj(Pq/p' q 2 — 1). More distant peaks (with lager value ui) 
have lager displacement, so that the interval occupied by 
the peaks contracts. At the same time, new peaks ap- 
pear within this interval, because with enlarging po the 
number of levels and the number of the possible tran- 
sitions between them increase. Distance between peaks 
decreases, and when it approaches H/t, the peaks begin 
to merge together. 

It is interesting to compare results of the study for the 
optical conductivity of ultrathin metal wires with anal- 
ogous results for ultrathin films. Divergences are asso- 
ciated with the different dimensionality of the systems. 
This is reflected in an essential difference in the energetic 
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Figure 7: Calculated frequency dependences of the real and 
imaginary parts of a a p (in e 2 /2aoh units) for Au wires of 
various diameter d. 



spectra and also in the fact that after calculation of qua- 
sicontinuous states, in the case of a wire, the summation 
over two numbers to, n remains, while in the case of a 
film, it remains over one number only (which numerates 
values of the x— component of the electron momentum). 
It is this fact that explains approximately one order lower 
height of the maximum in the frequency dependence of 
the conductivity of a wire compared to the case of a film. 
Indeed, 
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because 

(fcr c ) 



kl F <* 2nk° F /L. 



For L,p = 1 nm we obtain Re cr^ re /Re tr^f 10" 1 . 

As to the different position of the peaks, this may be 
completely explained by characteristic properties of spec- 
tra of the ID and 2D systems. 

The frequency dependences of Re a xx for the Al and 
Pb wires of diameter 1.6 nm are presented in Fig. 8. It 
is surprising that peaks in the conductivity of the Pb wire 
are absolutely absent. The reason of this is a small value 
of the relaxation time for Pb equal to r = 1.4 x 10~ 15 s, 
so that width of the peaks H/t = 0.44 eV. In this respect, 
Al, with fi/r = 0.08 eV, holds an intermediate position 
between Au and Pb. For calculations we use values of 
the relaxation time for bulk metals taken from [3|| . 

In spite of an absence of peaks in the frequency depen- 
dence of the conductivity for the Pb wire, its maximum 
may be found in such a way as the position of the max- 
imum height peak in the conductivity of the Au wire 
was determined above, with the difference that this time 
?ip = 4: 



(°1,4 



2.1 eV. 



2m e pl 

This value agrees well with Fig. 8 taking into account 
the large width of the peaks. 

Surprisingly, the difference in the obtained results for 
Al and Pb (due to the different values of r) indicate size- 
frequency dependences, which can be expected for films 
and wires, inhomogeneous in thickness considered, for 
example, in Refs. |40l - l44l ]. If fluctuations of sizes in 1D- 
and 2D-systems lead to a strong effective reduction of 
t, experimental size dependences of conductivity are no- 
ticeably smoothed irrespective of a metal kind. 

We will devote select publication to the theory of trans- 
port in films and wires with rough surface. 

VI. CONCLUSIONS 



ductivity tensor for a quasi-homogeneous ultrathin metal 
film and wire are calculated within the particle-in-a-box 
model on the assumption that the component of the in- 
duced current with the wave vector equal to the wave 
vector of the electromagnetic field is dominating. 

Over infrared region the condition gL, qpo <C 1 is sat- 
isfied allowing us to express components of the conduc- 
tivity tensor in terms of the according small value. All 
non-diagonal components of the conductivity tensor are 
equal to zero in zero order of the expansion. They ap- 
pear in linear approximation. The important fact that 
the real part of the diagonal components is non-negative 
over all frequency range, with the guarantee Q > for 
the dissipation of energy, is proved. 

As a result of comparing the according components of 
the conductivity tensor for a film and a wire of the same 
thickness of order 1 nm, one order smaller value for a 
wire is obtained. In such a manner different density of 
states near the Fermi level manifests itself (it is greater 
for a film). It is found that the discrepancy between 
our results and the theory [22j increases and becomes 
substantial, as the characteristic small dimension of the 
system decreases. This discrepancy is associated with 
the strong dependence of the Fermi level on this dimen- 
sion for small values of order of the Fermi wavelength. 
This size dependence of the Fermi level has an "oscil- 
latory" form. Transmittance is calculated for a simple, 
well-defined model without fitting parameters. 
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Appendix A: The matrix elements 

The expressions for momentum projections in cylindri- 
cal coordinates have a form 



Pz = -ioh—, 
dz 



. d sin ip d 
Px = -ion { cos<£- 



dp p dip J ' 



(Al) 



p y = —ioh s sin <p 



d cos ip d 



dp p dip 



The conductivity tensor is introduced for the low- 
dimensional electron systems. Components of the con- 



Using (|23J) - (J2HD and ([AT]) , we have 
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Then using relation I' m {x) — ±ml m (x) / x^f I m ±i(x) from In a similar way we find 
[13, we obtain Eqs. ([37 ]) - (|55 ]) . 
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